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1. INTRODUCTION 

Nowadays, multilevel inverter (MLI) has been applied in medium-voltage and high-power 
applications such as active filters, electrical motor drives and photovoltaic grid-connected system [1 ]-[7]. 
Multilevel inverter has several advantages such as lower switching losses, lower voltage stress on power 
switches, and lower electromagnetic interference (EMI) compared to two-level high-frequency pulse-width 
modulation (PWM) inverter [8], [9]. There are three main MLI topologies which are diode-clamped, flying 
capacitor and cascaded H-bridge [6]. Among them, cascaded H-bridge multilevel inverter (CHBMLI), which 
has the benefits of flexibility and modularity, has gained increasing attention in wide range of applications 
[6], [8]. CHBMLI is controlled by applying a set of optimum switching angles to synthesize a near sinusoidal 
staircase output voltage waveform. However, the optimum switching angles that can provide an output 
voltage waveform with low total harmonic distortion (THD) are not easy to be determined [8], [9]. 

Several methods to control the MLI have been reported in the past. Sinusoidal PWM and space 
vector PWM are the common high-frequency PWM methods to control the MLI [6], [10]. However, high 
switching loss is the main drawback in both methods. Another method is the fundamental-frequency PWM 
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method which has lower switching losses [6]. Optimal minimization of total harmonic distortion (OMTHD) 
is one of the fundamental-frequency PWM method that can minimize the THD [11]. However, it cannot 
guarantee that the low-order harmonics are minimized. Selective harmonic minimization pulse-width 
modulation (SHMPWM) is the common fundamental-frequaency PWM method that can be employed to 
obtain the optimum switching angles without having the problem of high switching losses [12]. SHMPWM 
can minimize the low-order harmonics and maintain the fundamental component at the same time. In this 
method, a set of non-linear transcendental equations of the selected harmonics that consist of trigonometric 
terms is to be solved. Newton Raphson (NR) technique is a common iterative mathematical method used to 
solve the equations [13]. However, NR is highly dependent on good initial switching-angle guesses and it 
often provides solution for a certain range of modulation index only. Another technique to solve the non- 
linear equations is by converting the non-linear equations to polynomial equations and then solving them 
using resultant theory (RT) [14], [16]. While this technique can provide a wide range of solutions, it is 
computationally complex and time-consuming for high-dimension switching-angle calculation. Hence, NR 
and RT are rather difficult to be implemented. Another method to optimize the MLI switching angles is by 
employing soft-computing approach [17]. The advantage of soft-computing approach is that the algorithm 
usually does not require a good guess of initial switching angles. Genetic algorithm (GA), which is a well- 
known soft-computing method, has been applied successfully in a wide range of applications, e.g. network 
routing and image processing [18]. It has also been employed to optimize the switching angles of MLI [19]. 
However, GA is easily trapped in the local optima due to the absence of exploration and exploitation 
abilities. Thus, GA has lower probability to find the global minima [20]. 

Recently, a nature-inspired soft-computing algorithm known as grasshopper optimization algorithm 
(GOA) has been proposed to solve optimization problem [21]. In last few years, GOA has been applied in 
several applications such as electrical characterization of proton exchange membrane fuel cells stack, power 
quality enhancement in an isolated microgrid, color image multilevel thresholding, meta-matching approach 
for ontology alignment, voltage control of switched reluctance generator and optimal reconfiguration of PV 
array [20], [22]-[26]. The application of GOA in THD minimization of power electronic converters has rarely 
been reported. Therefore, GOA is proposed in this paper to determine the optimum switching angles for the 
MLI. Unlike NR, GOA could find optimum solutions without requiring a good initial switching angle guess. 
The GOA based SHMPWM has been implemented and analyzed using MATLAB, whilst the switching-angle 
solution is compared to that obtained using NR and GA method. A Simulink/PSIM co-simulation model has 
also been developed to evaluate and verify the dynamic performance of a 3-phase 9-level CHBMLI with 
GOA optimized switching angles under dynamic modulation index demand. 


2. IMPLEMENTATION OF GOA-SHMPWM FOR 3-PHASE 9-LEVEL CHBMLI 

A CHBMLI is constructed by a series of H-bridge circuits. Each H-bridge circuit consists of a DC 
voltage source and four power semiconductor switches that are able to produce voltage levels +Va, O or -Vac. 
A (2k+1)-level of staircase output phase voltage waveform can be produced by a combination of k number of 
H-bridge circuits. Figure l(a) shows the construction of one of the phases and Figure 1(b) shows the output 
phase voltage waveform synthesized from four H-bridge circuits. In this paper, the phases of the CHBMLI 
are connected to form a balanced 3-phase Y-connection circuit with 120° of phase shift. Based on 
Figure 1(b), the output phase voltage waveform could be represented in Fourier series as (1) [7]: 


dc | cos(nay )+cos( naz )+cos(nar3 )+cos(nar4 ) sin(nat) 
m 





© AV, 
Von = X 
n 
n=l 


(1) 


where n is always an odd number and it also represents the n-th harmonics, Vac represents the magnitude of 
DC voltage source in each H-bridge circuit, and œ is the angular frequency of the fundamental harmonic. The 
switching angles a), @2, a3 and @4 are in unit radian and must satisfy the condition of 
0 < &ı < æ2 < &æ3 < q4 < 7/2 [19]. From (1), the n-th harmonic can be expressed as (2) [7]: 


4Vdc [ cos(nay )+cos(na2 )+cos(na3 )+cos(na4 )| 


nr (2) 





Vn = 


The objective of the SHMPWM is to minimize the undesired low-order harmonics and maintain the 
desired fundamental component of output voltage waveform. In a balanced 3-phase 9-level CHBMLI, 5th, 
7th and 11th harmonics are chosen to be minimized while the first harmonic is maintained at the desired 
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value. The triplen odd harmonis are not required to be minimized because they are eliminated naturally in the 
3-phase system. In this paper, GOA is applied in SHMPWM to optimize the switching angles. 


H-Bridge Circuit 





Figure 1. Block diagram showing a 3-phase 9-level CHBMLI and typical phase voltage waveform 


GOA is a population-based nature-inspired algorithm which is based on the swarm behavior of 
grasshoppers. The algorithm mimics the movements of the grasshoppers in migration and food foraging 
process. The mathematical model for implementing the movement of grasshoppers is given by (3) [21]: 





d d 
N XX: 
xf =c D eE f ed) j ʻi ve 
j=l dij 
j#i (3) 


where ubd is the upper bound in the d-th dimension, Ibd is the lower bound in the d-th dimension, Ta is the 
best target value and dj; is the distance between i-th grasshopper and j-th grasshopper. The c is a decreasing 
coefficient which is expressed as (4) [21]: 


CSC (iter Cae oan) / iter in) 


(4) 


where Cmax 1s the maximum value of the coefficient, Cmin 1s the minimum value of coefficient, iter is the 
current number of iteration and itermaxis the maximum number of iterations. The function s is the social force 
used to decide the movement of grasshopper and it is presented as (5) [21]: 


s= fe -e (5) 


where r is the normalized distance between i-th and j-th grasshoppers, / is the attractive length scale and f is 
the intensity of attraction. An objective function (OF), which is used in the GOA-SHMPWM to minimize the 
undesired harmonics and maintain the desired fundamental component, is adapted from [19] as (6): 


V j 1| 50 1 {50 : 1/50 
Open. 1) eel2 >) 222) 42) 


where Vp = (4kKMVj,)/z is desired fundamental harmonic that is controlled by the modulation index, M, and 
Vi, Vs, Vz and V7; are fundamental, 5th, 7th and 11th harmonics of the phase voltage waveform, respectively. 
The first term of (6) is used to regulate the desired fundamental harmonics, whilst the same time the second, 
third and forth terms are used to reduce the Sth, 7th and 11th harmonics, respectively. The objective of the 
first term is to limit the relative error between the Vp and V; by 1%. For the second, third and forth terms, the 
Sth, 7th and 11th harmonics are kept under 2% of the fundamental harmonic. Therefore, all desired 
conditions could be controlled with the proposed OF while the optimum solution for all modulation index 
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could be determined by using the implementation of GOA. Based on (2), Vi, Vs, Vz and V;; can be expressed 
as (7): 


= “ade 
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GOA-SHMPWM has been implemented using MATLAB to calculate the optimum switching 
angles, while the implementation of GODA-SHMPWM is illustrated in Figure 2. In the GOA-SHMPWM 
implementation, the parameters used for GOA are as follows: / = 1.5, f= 0.5, Cmin = 0.00001, Cmax = 0.5, N = 
100 and itefrmax = 100. Switching-angle calculation is performed for a 3-phase 9-level CHBMLI for 
modulation index range from 0.01 to 1.00 in step size of 0.01. 








Initialize population Xi (i = Calculate the OF value for each 


START 1,2, ...,n),Cmax, Cmin, l, f and M = 0.01 search agent using Eq. (7). 


itermax Ta = the best search agent. 





Normalize the distance between search agents in [1,4]. 


Update the position of each search agent using Eq. (4). 
Bring the search agent back if it goes outside the boundaries. 


iter = iter +1 
Update Ta if there is a better solution. 


Y 
Td = Global Minimum <a > 2 


Figure 2. Flowchart of GOA-SHMPWM implementation in a wide modulation index range 









In order to validate whether the CHBMLI controlled using GOA-SHMPWM optimum switching 
angles is capable to respond to the change of modulation index demand, a dynamic co-simulation of 
MATLAB/SIMULINK-PSIM is implemented on 3-phase 9-level CHBMLI. Figure 3 shows the Simulink 
model of switching-angle generator for 3-phase 9-level CHBMLI, whilst Figure 4 shows the block diagram 
of 3-phase 9-level CHBMLI simulation model and its PSIM implementation. Prior to the simulation, the 
switching angles that are optimized by GOA-SHMPWWM are stored in look-up tables as shown in Figure 3. 
During simulation, the optimized switching angles are interpolated from the look-up tables based on the 
dynamically changing modulation index demand. The switching angles are distributed into 24 signals (g1 to 
g24) through gate signal generator. By using the SimCoupler interface as shown in Figure 3, the 24 signals 
are linked to the PSIM 3-phase 9-level CHBMLI simulation model as shown in Figure 4. Each H-bridge 
module employs a 12V DC voltage source and the CHBMLI is switched at 50Hz. Output voltage waveforms 
are generated from the simulation model based on the modulation index demand. 
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Figure 3. Simulink model of switching-angle generator for 3-phase 9-level CHBMLI 
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Figure 4. Block diagram of 3-phase 9-level CHBMLI simulation model and its PSIM implementation 


3. RESULTS AND DISCUSSION 

Switching angles obtained from GOA are benchmarked against those obtained from NR technique 
and GA, as shown in Figure 5. Compared to NR as shown in Figure 5(a), GOA in Figure 5(b) has a wider 
modulation index range of optimum switching-angle solutions. Figure 6(a) shows the minimum OF achieved 
by GA and GOA for the modulation index range investigated. As shown in Figure 6(a), GOA mostly 
achieves lower OF compared to GA in a wide range of modulation index. As shown in Figure 6(b), 38% of 
the modulation index range achieves a minimum OF of 107 and below by using GA, and 38% of the 
modulation index range achieves a minimum OF of 10° and below by using GOA. This result suggests that 
GOA has a higher probability than GA to reach global minima in the CHBMLI optimization search space 
and the undesired 5th, 7th and 11th harmonics are nearly eliminated. As a result, GOA is able to produce 
higher accuracy of optimum switching angles than GA. Figure 6(c) shows the phase voltage THD and line- 
to-line voltage THD of the CHBMLI that achieved by GOA. The line-to-line voltage THD is always lower 
than the phase voltage THD because of the natural elimination of triplen harmonics in line-to-line voltage 
waveform. Based on Figure 6(c), the lowest phase voltage THD that can be achieved by GOA is 9.65% for 
M = 0.82. 
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Figure 5. Switching angles derived using (a) NR, (b) GA, (c) GOA techniques 
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Figure 6. MATLAB anslysis: (a) OF (GOA vs GA), (b) CDF (GOA vs GA), (c) Voltage THD (GOA) 


Table 1 shows the selected modulation indexes demand for the dynamic co-simulation, whilst 
Figure 7 shows the stepped modulation index demand, phase voltage waveforms and line-to-line voltage 
waveforms obtained from MATLAB/Simulink-PSIM dynamic co-simulation. The modulation index is varied 
at every 0.ls time interval. From 0 to 0.4s, when the modulation index is increased from 0.13 to 0.82, the 
voltage level increases from 3-level to 9-level and the peak phase voltage increases from 12V to 48V. At 
0.4s, the modulation index is decreased from 0.82 to 0.27, the voltage level decreases 9-level to 5-level and 
the peak phase voltage decreases from 48V to 24V. Due to the three-phase configuration, a 120° phase shift 
with respect to the adjacent leg is observed in the output phase voltage and line-to-line voltage waveforms. 
These results verify that the CHBMLI controlled using the GOA-SHMPWM optimized switching angles is 
capable to respond the change of modulation index demand. 

Close-up view of phase and line-to-line voltage waveforms for each modulation index demand are 
shown in Table 2. Fast Fourier Transform (FFT) are applied in the phase and line-to-line voltage waveforms 
for each modulation index to show the harmonic contents. Table 2 shows the voltage waveforms and FFT 
analysis for M = 0.13, M = 0.27, M = 0.41 and M = 0.82, respectively. The magnitudes of V; are 7.99V, 
16.54V, 25.10V and 50.10V, which are very close to Vp of 7.95V, 16.50V, 25.06V and 50.11¥V as listed in 
Table 1, respectively, for M = 0.13, M =0.27, M= 0.41 and M = 0.82, respectively. In addition, the 
magnitude of fundamental harmonic of the line-to-line voltage waveform for each modulation index is 
approximately V3xV;. For M = 0.27, M =0.41 and M = 0.82, the chosen 5th, 7th and 11th harmonics in the 
phase and line voltage waveforms are nearly eliminated. This is reasonable since the OF related with these 
modulation indexes are successfully minimized, which can be confirmed in Figure 6(a). Neverthelessly, the 
Sth, 7th and 11th harmonics of the voltage waveforms for M = 0.13 are only reduced as much as possible due 
to the low modulation index of output voltage waveform that has the lowest number of voltage levels. 
Table 3 compares the MATLAB and dynamic simulation results of fundamental harmonic, phase voltage 
THD and line-to-line voltage THD for each modulation index demand. It can be observed that the results 
between MATLAB and dynamic simulation are very close. 
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In order to demonstrate that the CHBMLI controlled using GOA-SHMPWM switching angles is 
able to produce a steady output waveform during the load change, the load impedence is dynamically varied 
in the simulation. Figure 8 shows the change of phase voltage waveform and load current waveform during 
the dynamic load change at 0.1s for modulation index M = 0.82. Before 0.1s, the impedence is 1Q and power 
factor is 0.81 while the peak value of the phase voltage and load current is 48V and 51.91A, respectively. 
After 0.1s, the impedence is 0.72Q and power factor is 0.57, whilst the peak phase voltage and load current is 
48V and 74.86A, respectively. This figure shows that the peak value of the load current increases during the 
load change and the phase voltage waveform still remain the same after the load change. 


Table 1. Optimum switching angles and number of phase voltage levels of selected modulation index used 
for 3-phase 9-level CHBMLI MATLAB/Simulink-PSIM dynamic co-simulation 


M Vp (V) a, (deg.) Qn (deg.) a3 (deg.) a4 (deg.) Num. of voltage levels 

0.13 7.95 58.46 90.00 90.00 90.00 3 

0.27 16.50 41.89 67.77 90.00 90.00 5 

0.41 25.06 25.05 51.62 64.31 90.00 7 

0.82 50.11 8.63 19.22 34.69 58.34 9 
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Figure 7. Stepped modulation index demand at every 0.1s interval, phase voltage waveforms and line-to-line 
voltage waveforms obtained from MATLAB/Simulink-PSIM dynamic co-simulation 
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Figure 8. Dynamic load change is applied at 0.1s (IZI is load impedance; M = 0.82) 
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Table 2. Closed-up view of simulated output voltage waveforms and FFT analysis 


M Dynamic Co-Simulation Results 
Phase and Line Voltage Waveforms Phase Voltage FFT Analysis Line Voltage FFT Analysis 
[voltage (V) vs time (s)] [harmonic amplitude (V) vs freq. (Hz)] [harmonic amplitude (V) vs freq. (Hz)] 
0.13 100 100 





Vi Vs Vo  V1(7.99V) 



















































vii | 
50 50 foe eee eee Cee] eee ee — B E me eee See (oer oat ee 
25 irt cceceEs 25 
0 0.02 0.04 0 250 500 250 1k 0 
0.27 100 
75 
0.41 Se 
O n O Ea 
0.82 100 





V1 (50.10V) 
Vu l 

















Table 3. Fundamental harmonic, phase voltage THD and line-to-line voltage THD comparison between 
MATLAB calculation result and PSIM simulation result 


_ Comparison CW MATLAB Analysis Result Dyn. Co-Sim. Result CD fference __ 

V; (V) 7.99 7.99 0.00 

M =0.13 Vn THD (%) 76.13 76.11 0.02 
Vine THD (%) 32.41 32.38 0.03 

V; (V) 16.54 16.54 0.00 

M =0.27 Vn THD (%) 52.24 52.25 0.01 
Vine THD (%) 19.74 19.74 0.00 

V; (V) 25.11 25.10 0.01 

M=0.41 Vn THD (%) 45.97 45.97 0.00 
Vine THD (%) 13.74 13.75 0.01 

V; (V) 50.11 50.10 0.01 

M = 0.82 Vn THD (%) 9.65 9.66 0.01 
Vine THD (%) 5.80 5.80 0.00 


4. CONCLUSION 

In this paper, the performance of GOA for optimizing switching angles applied to dynamic 
simulation of 3-phase 9-level CHBMLI have been presented. The switching angle results show that GOA has 
a wider range of optimized switching-angle solutions than NR technique. The OF and CDF results show that 
GOA is able to determine the optimum switching angles with higher accuracy compared to GA. Thus, GOA 
has superior performance over NR and GA in the MLI control. In addition, the MATLAB/Simulink-PSIM 
dynamic co-simulation result verifies that the CHBMLI can respond to the change of modulation index 
demand and generate output voltage waveforms with minimized selected harmonics and desired fundamental 
component. Furthermore, the simulation result of dynamic load change validates that CHBMLI with GOA- 
SHMPWM switching angles is able to produce a steady output phase voltage under dynamic load change. 
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